Momentum signatures for Schwinger pair production in short laser pulses with a 

sub-cycle structure 
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We investigate electron-positron pair production from vacuum for short laser pulses with sub-cycle 
structure, in the nonperturbative regime (Schwinger pair production). We use the non-equilibrium 
quantum kinetic approach, and show that the momentum spectrum of the created electron-positron 
pairs is extremely sensitive to the sub-cycle dynamics - depending on the laser frequency oj, the pulse 
length r, and the carrier phase <f> - and shows several distinctive new signatures. This observation 
could help not only in the design of laser pulses to optimize the experimental signature of Schwinger 
pair production, but also ultimately lead to new probes of light pulses at extremely short time scales. 

PACS numbers: 12.20.Ds, ll.15.Tk, 42.50.Xa 



Electron-positron pair production due to the instabil- 
ity of the quantum electrodynamics (QED) vacuum in an 
external electric field is a remarkable non-perturbative 
prediction of QED [l|, 0, 0] that has not yet been directly 
observed. Significant recent advances in laser technology 
have raised hopes that the required critical field strength 
of ~ 10 16 V/cm may soon be within experimental 
reach [1, [H, [f|, either in optical high-intensity laser fa- 
cilities such as Vulcan or ELI [3] or in X-ray free elec- 
tron laser (XFEL) facilities . Observation of this elu- 
sive phenomenon in the non-perturbative domain would 
complement the perturbative multi-photon pair produc- 
tion seen at SLAC E144 using nonlinear Compton scat- 
tering [jj. Moreover, it would represent a significant 
advance in our understanding of non-perturbative phe- 
nomena in quantum field theory, with potentially impor- 
tant lessons for related phenomena such as Unruh and 
Hawking radiation. The Schwinger mechanism has also 
been used to study various non-perturbative phenomena: 
e.g., string-breaking in the strong interactions [13], pair 
production in supercritical fields [11 1, n eutrino produc- 
tion in a fermionic density gradient 12J, and saturation 
in heavy-ion collisions [13]. Since the basic physics is 
quantum tunneling, the effect is exponentially weak, and 
so it is important to search for distinctive signs that 
might facilitate its detection. Here we consider a real- 
istic laser pulse with sub-cycle structure, and find dis- 
tinctive new signatures in the momentum distribution of 
the produced pairs. We also explain these signatures by 
rela ting the non-equilibrium quantum kinetic approach 
[Ti . Il5l. fl6l. UtI [Til , flip to the quantum mechanical scat- 
tering description [53, El, H3] ■ 

The original estimates [l], 0, 0] assumed a constant 
and uniform external electric field, but realistic ultra- 
strong fields are realized in short pulse, focussed lasers. 
We concentrate here on the time dependence of the elec- 
tromagnetic field and neglect spatial variations, assum- 
ing that the spatial focussing scale is much larger than 
the Compton wavelength. This approximates the exper- 
imental situation of two counter-propagating short laser 



pulses, generating a standing- wave electric field which is 
approximately spatially homogeneous in the interaction 
region, such that E(t) = (0, 0, E{t)) with (see Fig. [J) 



E(t) 



E cos(cj£ + (f>) exp ( — 



(1) 



Here lo is the laser frequency, r defines the total pulse 
length, and <fi is the 'carrier phase' (carrier-envelope 
absolute phase). We are motivated to investigate the 
carrier-phase dependence of the Schwinger mechanism 
by the sensitive carrier-phase dependence of strong-field 
ionization experiments in atomic, molecular and optical 
(AMO) physics [23]. It is convenient to introduce the pa- 
rameter a = ujt as a measure of the number of oscillation 
cycles within the Gaussian envelope pulse. This type of 
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FIG. 1: Shape of the electric field Eq. ([T]), for carrier-phase 
<j> — 0, when passing from a = 3 (dotted line) to a = 4 
(dashed line) to a = 5 (solid line). The pure Gaussian field 
(dashed-dotted line) is given as reference. 

electromagnetic field configuration can be represented by 
a time-dependent vector potential A(t) = (0,0, A{t))\ 



A{t) 



t . a \ 



+c.c. (2) 



Since Schwinger pair production is a time-dependent 
nonequilibrium process, quantum kinetic theory provides 



2 



an appropriate framework [17J, [lj, |2J, |25|, |2fj, |27J, |2& |2jJ, 
l3C| . (Below, we relate this to the widely-used WKB ap- 
proach). Within the quantum kinetic approach, the key 
quantity is the single-particle momentum distribution 
function f(k,t), which satisfies a non-Markovian quan- 
tum Vlasov equation including a source term for electron- 
positron pair production. We stress that f(k,t) is physi- 
cally meaningful as the distribution function of real parti- 
cles only at asymptotic times t — > ±00, when the electric 
field vanishes. As demonstrated in [24|, the field-current 
feedback due to Maxwell's equation can be neglected in 
the subcritical field strength regime, Eo < 0.1-E C r- Thus, 
we need to solve just one integro-differential equation: 



d/(M) = 1 eE(t)e ± 
dt ' 2 ^(k,t) 



1 M e -m^ 

-00 L0 2 (k,t') 



i-2/(M') 









x cos 
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c,r) 









(3) 



Here e is the electric charge; k = (k±,k^) is the canon- 
ical three-momentum; the kinetic momentum along the 
field is defined as p\\(t) = fc|| — eA(t); e\ = m 2 + k 2 is 
the transverse energy squared, and u) 2 {k,t) = e\ + p 2 (t) 
characterizes the total energy squared. Note that Eq. J3]) 
is valid for Dirac particles (QED), however, there exists 
a very similar equation for scalar particles (sQED) which 
takes the statistics into account as well [171 llSl [l9j . 

The appearance of, and interplay between, a total of 
four scales - electron mass to, applied electric field Eq, 
laser frequency ui, and total pulse length r - makes 
the electron-positron pair production process in such a 
laser field rather complicated [3l|, and suggests that the 
physics will not simply depend on the Keldysh parame- 
ter 7 = muj/eEo. Schwinger pair production in a pulsed 
electric field Eq. |TJ), with zero carrier phase (0 = 0), has 
been studied using WKB j32|. In the nonperturbative 
regime, 7 < 1, the momentum spectrum is: 



d 3 V 
dk* 



CXp — 7T 



E C i- 



1 

^E~o 



YK + kl (4) 



with 7 2 = (1 + l/cr 2 )7 2 . We will show that several inter- 
esting properties of the momentum distribution function 
for short pulses with many cycles per pulse are not cap- 
tured by WKB. 

We numerically integrate the quantum Vlasov equation 
Eq. ([3]) to obtain the asymptotic distribution function 
f(k,oo) = d 3 7 3 /dfc 3 , for various values of the laser fre- 
quency uj, choosing a temporal width r = 2 ■ 10 _4 cV _1 , 
which corresponds to a total pulse length of several times 
1CP 19 s. This lies in the anticipated range of experimen- 
tal parameters of future XFELs or may become realizable 
with higher harmonics or secondary-beam generation of 
optical lasers. We concentrate on the dependence on the 
longitudinal momentum fc|| , setting k± = 0. For the more 
realistic case of additional spatial pulse inhomogeneities, 
we expect also nontrivial information encoded in all spa- 
tial momentum directions. 



The most dramatic new effect is that the momentum 
distribution shows distinctive oscillations, with an oscil- 
lation scale set by the laser frequency, as shown in Fig. [2j 
The distance between successive peaks in f(k, 00) is given 
by u. This oscillatory behavior becomes pronounced 
when a > 4, and the amplitude of the oscillations in 
the distribution function increases further as we increase 
the number of cycles within the pulse beyond 4. 
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FIG. 2: Asymptotic distribution function f(k, 00) for k± — 0, 
Eq = 0.lE cr and <t> = when passing from a = 3 (dotted line) 
to a — 4 (dashed line) to a = 5 (solid line). Note that for 
cr = 3 the canonical momentum value fey = corresponds to 
a kinetic momentum P||(oo) ~ 75keV, whereas for a = 5 the 
value fcy = is equivalent to p\\ (00) m 0. 

To understand the physical origin of these oscillations, 
we recall that for a spatially uniform but time dependent 
electric field, the pair production process can be viewed 
as a one dimensional quantum mechanical scattering 
problem, with 'potential' given by -uj 2 (k,t) [13, El, E3| . 
The number of produced pairs is related to the reflection 
coefficient for this over-the-barrier scattering problem. In 
fact, this scattering picture is completely equivalent to 
the quantum kinetic approach [H^, [H[, and so we can 
interpret these oscillations as being due to resonances 
in the scattering problem as the shape of the potential 
—<jj 2 {k,t) changes with the momentum k. These oscilla- 
tions are missed when employing the WKB approxima- 
tion to the scattering problem on which Eq. ([4]) is based, 
but can be seen clearly in an exact numerical integration 
of the scattering problem, or as we have shown here, in 
an exact numerical integration of the quantum Vlasov 
equation. This physical picture explains why the spac- 
ing between the peaks is the laser frequency, and also 
explains the sensitive dependence on the other shape pa- 
rameters, such as cr. 

In fact, there is an even more distinctive dependence 
on the carrier phase, (ft, upon which the form of the scat- 
tering potential —u> 2 (k, t) is extremely sensitive. The car- 
rier phase dependence is difficult to discuss in the WKB 
approach, because a nonzero carrier phase breaks the 
E(t) = E(—t) symmetry of the pulse shape, which in turn 
makes the imaginary time treatment of the WKB scatter- 
ing problem significantly more complicated [2l|. But in 
the quantum kinetic approach, the carrier phase causes 
no computational problems; it is just another parame- 
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FIG. 3: Asymptotic distribution function f(k, oo) for k± — 
for a — 5, Eo = 0.1E cr and <f> — — 7r/4. The center of the 
distribution is shifted to py (oo) w 102 keV. 
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FIG. 4: Asymptotic distribution function f(k, oo) for fcx = 
for a — 5, So = 0.1-E C r and <f> = —ty/2. The center of the 
distribution is shifted to P||(oo) ~ 137 keV. 



ter. We have found that the introduction of the carrier 
phase makes the oscillatory behavior in the longitudinal 
momentum distribution even more pronounced. This is 
shown in Figs. [3] and [H where the momentum distribu- 
tion function is plotted for <j> = — 7r/4 and 4> = — 7r/2. We 
see that for the same values of the other parameters, the 
oscillatory behavior becomes more distinct as the phase 
offset increases. The most distinctive momentum signa- 
ture, however, is found for 4> = ~ 7r /2, when the electric 
field is totally antisymmetric. In this case, the asymp- 
totic distribution function f(k,t) vanishes at the minima 
of the oscillations, as shown in Fig. |H This feature also 
has a direct analogue in the scattering picture: for an 
antisymmetric field, the gauge potential Eq. J2]) is sym- 
metric and so is the scattering potential well — LJ 2 {k,t). 
In this case, perfect transmission is possible for certain 
resonance momenta, corresponding to zero reflection and 
thus zero pair production. Also note that the center of 
the distribution shifts from py (oo) = to a non-zero value 
again. These carrier-phase effects provide distinctive sig- 
natures, strongly suggesting a new experimental strategy 
and probe in the search for Schwinger pair production. 

These momentum signatures can also be understood 
in a quantum-mechanical double-slit picture, which has 
first been developed in the context of above-threshold 
ionization with few-cycle laser pulses (35|: in this pic- 



ture, the oscillations are fringes in the momentum spec- 
trum that result from the interference of temporally sep- 
arated pair creation events. The fringes are large for 
4> = — 7r/2, since then the field strength has two peaks of 
equal size (though opposite sign) which act as two tempo- 
rally separated slits. Moving the carrier phase away from 
4> = —n/2 corresponds to gradually opening or closing 
the slits, resulting in a varying degree of which- way in- 
formation and thus a varying contrast of the interference 
fringes. A quantitative consequence of this double-slit 
picture is that the width of the envelope of the oscilla- 
tions in the distribution function is related to the tem- 
poral width of the slits. The width of the envelope of 
oscillations thus also becomes a probe of the sub-cycle 
structure of the laser. 

To complete the physical picture, we consider the over- 
all envelope of the longitudinal momentum distribution, 
again for <j> = 0, averaging over the rapid oscillations. 
When there are more than three cycles per pulse (a > 3), 
the peak of the momentum distribution is located near 
P||(oo) = 0, whereas for a < 3 the peak is shifted to a 
non-zero value. Furthermore, the Gaussian width of the 
employed WKB approximation Eq. ([4]), which scales with 
v / eSo/7, is obviously somewhat broader than the true 
distribution, as is shown in Fig. EJ We can quantify this 
discrepancy in the width, by extending the WKB result 
beyond the Gaussian approximation inherent in Eq. - 
We use the results from (3(|, where it has been shown 
that the WKB instanton action 5g in scalar QED can 
also be applied for spinor QED. Within this approach, 
d 3 V/dk 3 ~ exp(— 2iSV), where the instanton action can 
be defined in the complex t— plane as a contour integral: 

2Sj:^i^ yjm 2 + k\ + jjfej - eA{t)] 2 At , (5) 

with the path Y around the branch cut. After a change 
of variable, from t to T = —A(t)/Eo, we can expand 
the instanton action in powers of (m/eEor) 2 and 
(fey / eEgr) 2 . The result of this improved WKB expression 
to third order is plotted as the dashed line in Fig. [5j while 
the Gaussian WKB formula Eq. is shown as a dotted 
line. We see that exp(— 25^) fits very well the averaged 
envelope of the exact momentum distribution, while the 
Gaussian approximation Eq. (@]) is significantly broader. 
Neither WKB estimate sees the oscillatory structure of 
the momentum distribution. Furthermore, the oscilla- 
tory behaviour of the momentum distribution function 
is found also for scalar particles in the quantum kinetic 
framework of sQED, with the statistics playing a crucial 
role: the averaged envelope of the momentum distribu- 
tion functions is identical, however, at momentum values 
where QED predicts a local maximum in the momentum 
distribution, sQED predicts a local minimum, and vice 
versa. 

To conclude, we point out that the momentum distri- 
bution signatures described here, as well as their carrier- 
phase and oscillation-number dependence, provides for a 
new handle on the first detection of Schwinger pair pro- 
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FIG. 5: Comparison of the asymptotic distribution function 
f(k, oo) for k± — (oscillating solid line) with the prediction 
of Eq. ([4]) (dotted line) and the improved WKB approxima- 
tion based on an expansion of Eq. ([U (dashed line) for a = 5, 
E = O.lScr and = 0. 

duction and its quantitative exploration, once laser field 
strengths have become sufficiently strong. In practice, 
these momentum signatures can distinguish Schwinger- 
produced pairs from possible background events, e.g., in- 
duced by residual-gas effects. As electron spectrometers 
used in laser-plasma acceleration experiments can reach 



a resolution of better than 1%, the whole spectrum in- 
cluding the oscillations will be directly accessible. 

Analogously to strong-field ionization experiments (23l . 
EH, E3, S|, these new pair-production signatures may 
also serve as sensitive probes of sub-cycle structure in 
ultra-short laser pulses. In addition to the positron yield, 
our results suggest a number of new observables such as 
the peak position and the width of the momentum dis- 
tribution function and, most importantly, its potentially 
oscillatory structure. In particular, the characteristics 
of the oscillations can provide rather precise information 
about the carrier phase and the total pulse length. As the 
latter laser characteristics are difficult to control a priori 
in an absolute manner, these momentum signatures can 
serve as a tomograph of the laser pulse, providing for a 
unique means to verify and confirm design goals of future 
laser systems at highest intensities. 
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